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ABSTRACT 

Radio-loud active galactic nuclei are among the most powerful objects in the universe. In these objects, most of the emission comes 
from relativistic jets getting their power from the accretion of matter onto supermassive black holes. However, despite the number of 
studies, a jet’s acceleration to relativistic speeds is still poorly understood. 

It is widely known that jets contain relativistic particles that emit radiation through several physical processes, one of them being the 
inverse Compton scattering of photons coming from external sources. In the case of a plasma composed of electrons and positrons 
continuously heated by the turbulence, inverse Compton scattering can lead to relativistic bulk motions through the Compton rocket 
effect. We investigate this process and compute the resulting bulk Lorentz factor in the complex photon field of an AGN composed of 
several external photon sources. 

We consider various sources here: the accretion disk, the dusty torus, and the broad line region. We take their geometry and anisotropy 
carefully into account in order to numerically compute the bulk Lorentz factor of the jet at every altitude. 

The study, made for a broad range of parameters, shows interesting and unexpected behaviors of the bulk Lorentz factor, exhibiting 
acceleration and deceleration zones in the jet. We investigate the patterns of the bulk Lorentz factor along the jet depending on the 
source sizes and on the observation angle and we finally show that these patterns can induce variability in the AGN emission with 
timescales going from hours to months. 
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1. Introduction 

It is now widely known that AGN jets hold relativistic flows. 
The first evidence of this goes back to the 1960s with the inter¬ 
pretation of brightness temperatures of quasar radio emission ex¬ 
ceeding the Compton limit by Rees (1966). This was beautifully 
confirmed by the observation of superluminal motions with the 
achievement of the very long base interferometry (VLBI) tech¬ 
nique (Cohen et al. 1971). These superluminal events are only 
possible for actual speeds larger than 0.7c. Relativistic velocities 
are also required to avoid strong y - y absorption by pair pro¬ 
duction, the high-energy photons being able to escape thanks to 
beamed radiation (see Baring 1994). 

Relativistic flows are characterized by their bulk Lorentz fac¬ 
tor = (l -filj rather than their speed \4 with = V^/c. 
However, there are several pieces of evidence that this bulk 
Lorentz factor is not homogeneous throughout the flow. The spa¬ 
tial distribution of the relativistic motion in the jet is still a mat¬ 
ter of discussion. Two types of variations are possible: radial 
and longitudinal (or a combination of both). The variations of 
longitudinal bulk Lorentz factor have often been parametrized 
by power laws with an accelerating and/or a decelerating phase 
(Marscher 1980, Ghisellini et al. 1985, Georganopoulos & 
Marscher 1998, Li & Wang 2004, Boutelier et al. 2008). Even 
though an initial accelerating phase appears necessary to achieve 
relativistic speeds, decelerating flows have also been invoked at 
larger scales, for example, to unify BL Lacs and radiogalaxies. In 
the unification scheme, BL Lacs and FR I radio galaxies are the 
same type of objects seen at different angles. However, BL Lacs 
models with constant jet velocities need very high bulk Lorentz 


factors to produce the observed SED, which is in contradiction 
with the FR I models and the observed jet velocities at subparsec 
scales in the TeV BLacs Mrk 421 and Mrk 501 (Marscher 1999). 
In their work, Georganopoulos & Kazanas (2003) showed that 
deceleration of the jet allows photons emitted in the inner parts 
of it to be scattered by the upper parts. In this case, radiation 
from fast regions of the jet would be highly beamed and thus 
correspond to BL Lacs objects, while radiation from slower re¬ 
gions would be emitted in a wider cone and would correspond to 
radio galaxies. Implications of a bulk velocity structure for the 
observed spectral energy distribution (SED) has been studied by 
Yang et al. (2009). 

In addition, a radial distribution of velocity is possible and 
has been particularly studied in the case of a double-jet structure, 
the so-called spine/layer jet. Here too, the idea was proposed 
as a solution to the unification issue between BL Lac objects 
and radiogalaxies (see Chiaberge et al. 2000). In this framework, 
Ghisellini et al. (2005) explain the rapid variability of the TeV 
emission without requiring huge Doppler factors, and Tavecchio 
& Ghisellini (2014) are able to reproduce the SED of NGC 1275. 
Recent observations (Giovannini et al. (1999), Swain et al. 1998, 
Giroletti et al. 2004) also bring evidence of such structures. 

The idea of a two-flow structure was first proposed by Sol 
et al. (1989) for theoretical reasons. In this paradigm, the jet 
is assumed to be made of two components: a mildly relativistic 
sheath composed of e~ /p + and driven by magnetohydrodynam- 
ical (MHD) forces, which transports most of the kinetic energy, 
and an ultra-relativistic spine composed of e~/e + pairs, which is 
responsible for most of the emission . A detailed description of 
the formation of such a pair beam has been developed in Mar- 
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cowith et al. (1995) and following works (e.g., Renaud & Henri 
(1998), Henri & Sauge (2006), Boutelier et al. (2008)). In this 
model, the outer jet acts as an energy reservoir for the particles 
of the spine. Starting from an initial injection of some relativis¬ 
tic particles (possibly created in the surrounding of a rotating 
black hole), these particles will emit copious amounts of high- 
energy radiation, which will be converted into pairs. These pairs 
will be in turn continuously reaccelerated along the jet via the 
second-order Fermi process through the turbulence triggered by 
various instabilities in the outer MHD jet. Observations of dif¬ 
fuse X-ray emission in type 1 Fanaroff-Riley objetcs are in favor 
of this view of distributed particle acceleration rather than local¬ 
ized shocks (Hardcastle et al. 2007). 

In leptonic models, X-ray and gamma-ray emission is 
thought to arise from inverse Compton (IC) process on soft pho¬ 
tons. These photons can be provided by synchrotron emission 
(called the synchrotron self-Compton or SSC process), or by ex¬ 
ternal sources such as an accretion disk, the broad line region, 
or a dusty torus. All these sources will give a locally anisotropic 
photon field on the axis of the jet. Under these conditions, the 
emitted radiation will also be highly anisotropic, which produces 
a strong momentum transfer between the relativistic plasma and 
the emitted radiation, the so-called Compton rocket effect (see 
O’dell 1981, Melia & Konigl 1989 or more recently Renaud & 
Henri 1998). 

For a relativistic e~/e + pair plasma, this force is dominant 
and will drive the bulk motion to relativistic velocities. As is 
detailed in section 2, this effect will saturate when the velocity 
(or equivalently the bulk Lorentz factor) of the plasma reaches 
a characteristic value for which the net radiation flux vanishes 
in the comoving frame as a result of the aberration of the photon 
momentum. Thus, a plasma submitted to this radiation force will 
tend to reach this equilibrium velocity, which can be viewed as 
the average velocity of the photon “wind”. The Compton rocket 
effect can also be found under the name of Compton drag effect. 
Even though they are exactly the same mechanism, the drag de¬ 
nomination comes from an a priori assumption of a very high 
bulk Lorentz factor, higher than that imposed by the external ra¬ 
diation and the Compton rocket effect: in this case, the inverse 
Compton rocket (more appropriately called a “retrorocket” ef¬ 
fect in this case) will result in a deceleration of the flow or a 
limitation of its velocity (Sikora et al. (1996) and Ghisellini & 
Tavecchio (2010)). 

The Compton rocket effect has often been dismissed as a 
cause of relativistic motion because it is also a cooling pro¬ 
cess. This means that an isolated relativistic pair plasma will 
also be quickly cooled and will generally be unable to reach the 
high bulk Lorentz factors (« 10) needed to explain superluminal 
motions (Phinney 1982 and also Madau & Thompson (2000)). 
However, this objection is not valid in the two-flow model, since 
the relativistic pair plasma is supposed to be continuously re¬ 
heated by the surrounding MHD flow. 

Under these conditions, a pure electron-positron pair plasma 
can be coupled to the radiation field over a much larger distance. 
It will then stick to the equilibrium velocity (which is generally 
variable), until the radiation field weakens enough for the decou¬ 
pling to occur. The plasma will then essentially follow a ballistic 
motion at the terminal velocity, which depends on the location 
of the decoupling. In the following, we will only study the value 
of the equilibrium velocity, which depends only on the radiation 
field and not on the characteristics of the plasma in the Thomson 
regime. On the other hand, the location of the decoupling, and 
hence the terminal velocity, depends on these plasma character¬ 
istics, as is discussed in section 6. 


In summary, we investigate this paradigm and propose to 
study the evolution of the resulting bulk Lorentz factor (pre¬ 
sented in section 2) due to the Compton rocket effect in a com¬ 
plex photon field including three main external sources of soft 
photons present in quasars, the accretion disk, the dusty torus, 
and the broad line region (BLR) (presented in section 3). By 
computing accurately the equilibrium bulk Lorentz factor along 
the jet in the Thomson regime (sections 4 and 5), we can study 
what effect it might have on the observed emission (sections 5 
and 7). 

In the rest of the paper, the “jet” refers to the inner spine 
of the two-flow model which is subject to the Compton rocket 
effect and is the flow moving at relativistic bulk speeds. Primed 
quantities are expressed in the comoving frame and unprimed 
quantities are expressed in the lab frame (i.e., the external source 
frame of the AGN). 

2. T bu i k and equilibrium 

We assume a static bulk of relativistic leptons following an 
isotropic distribution in an anisotropic photon field. Owing to 
the Doppler effect, particles moving towards the main light 
source will scatter photons of higher energy and with a higher 
rate than those moving outwards. This will naturally lead to an 
anisotropic inverse Compton emission, most of it going back to 
the main photon source. This anisotropic emission will result in 
a transfer of momentum on the emitting plasma, the so-called 
Compton rocket effect first described by O’dell (1981). A hot 
plasma could be driven to relativistic bulk motion through this 
mechanism. The force depends on the anisotropy of the soft 
photons seen in the comoving frame. When relativistic motion 
is taken into account, the photon field in the comoving frame 
will be affected by the bulk Doppler factor, resulting in a more 
isotropic photon distribution, until the Compton rocket force 
vanishes. The plasma then reaches an equilibrium velocity, or 
equivalently, an equilibrium bulk Lorentz factor Y eq . 

In the Thomson regime, T eq depends only on the external 
radiation field through the Eddington parameters (/, H, K) 

J = L J'I Vs (£l s )dQ s dv s 

H = L f Iy s (Q s )fi s dQ s dv s (1) 

K =^f W^dn.dvs 

with I Vs the specific intensity of the emitting source, the 
solid angle, and 6 S = arccos p s the angle under which the source 
is seen by the pair plasma (see figure 1). 

In the Thomson regime, the saturation of the Compton rocket 
effects happens when the second Eddington parameter, //', van¬ 
ishes in the comoving frame (see Marcowith et al. 1997 for more 

J K 

details). With the factor £ = —-, one obtains the equilibrium 

2 H 

equation 

H' = H Y 2 eq (p 2 eq - 2£/3 eq + l) = 0, (2) 

whose solution is 

0 eq = (- (3) 
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It is interesting to note again that, as long as the plasma is 
hot, this result does not depend on the jet model in the Thomson 
regime, but only on the external photon field (see equation 
2). The Compton rocket effect will also take place in the 
Klein-Nishina (KN) regime, but in this case the computation of 
the equilibrium velocity is a bit more complex and depends on 
the energy distribution of the pair (see Renaud & Henri (1998)). 
Moreover, one can expect that the resulting T eq will not vary 
much from that computed in the Thomson regime as long as we 
are not deep in the KN regime, as the recoil in the Thomson 
regime is much more efficient. Then Y eq could be sensitive to 
KN correction for the most extreme objects (> lOOGeV; see 
Sect. 7). In the following, we suppose that conditions always 
meet the Thomson regime. The goal of this paper is to compute 
the resulting equilibrium bulk Lorentz factor T eq , in a complex 
environnement, taking into account the angular and spectral 
distribution of various sources of photons in a realistic model of 
AGNs. 


3. Modeling the AGN 

We will consider the effect of several possible external sources 
of photons, namely the accretion disk, the dusty torus, and the 
broad fine region (see figure 1). We note that because the syn¬ 
chrotron radiation is produced in the comoving frame with a zero 
net flux, it does not interfere with the bulk motion as long as 
the SSC is treated in a local approximation. This could change, 
however, if the particles scatter synchrotron photons produced 
in other parts of the jet, but this problem is much more complex 
since it involves the knowledge of the whole structure of the jet. 
We will not address this issue in this work. 


external emission towards the observer’s fine of sight. Since the 
object is seen under a certain angle, the axisymmetry is always 
broken with respect to the line of sight. 

Each slice of the discretization is then described by four 
numbers: 


- Us = cos 0 S with 6 S the angle between the incoming light 
wave and the jet axis (see figure 1), 

- d£2 the solid angle under which it is seen from the altitude Z, 

- its temperature T , 

- its emissivity s. 

The slice emission characterized by these numbers is given 
by its specific intensity I v (erg s~ l cm ~ 2 sr~ x Hz ~defined as 
the emitted energy d E by normal surface dX, time c it, frequency 
band dv, and solid angle d Q: dE = I v dLdtd£ldv (Rybicki & 
Lightman 1979). The specific intensity of a graybody is given 
by Planck’s law (equation 4): 


I v = s- 


2hv 3 


exp 




( 4 ) 


A blackbody is a graybody with an emissivity s = 1. 

Here below, we detail the computation of each photon 
source. 


3.2. Standard accretion disk 



4 





BLR . 

N # * ’ 


/ 


, 

accretion : black hole 
disk 




Fig. 1: The big picture: edge-on view of the global model geom¬ 
etry (not to scale) with the accretion disk, the dusty torus, and 
the broad line region (BLR). i 0 b s is the observer’s viewing angle 
and 6 S is the angle between the incoming radiation from a source 
and the jet axis. 


3.1. Discretization of the sources 

The anisotropy of the photon sources will be taken into account 
in the numerical scheme by slicing the different sources into 
a set of small independent parts modeled as graybodies in 
radiative equilibrium, i.e., with a Planck spectrum but with a 
possible smaller emissivity. This discretization is done in three 
dimensions. Even for an axisymmetric source, an azimuthal dis¬ 
cretization is still required to accurately compute the Compton 


The accretion disk is assumed to be an optically thick standard 
accretion disk as described by Shakura & Sunyaev (1973) ex¬ 
tending from R in to R out . Each point of the disk is a blackbody 
with a temperature given by the distance to the central black hole 
(BH), assumed to be non-rotating, following the relation 


Tdiskir) - 


3 GMM 1 ( 

8 7rcr r 3 y 



( 5 ) 


with G the gravitational constant, cr the Stefan-Boltzmann 

GM 

constant, M the BH mass, M the accretion rate, Rs = 2—— the 

c 2 

Schwarzschild radius, and r the distance from the center of the 
BH. 


To model this accretion disk, we sliced it in different parts, 
each being a blackbody with its own temperature T disk (r, cp). The 
discretization follows a logarithmic scale along r and a linear 
scale along p. Given an altitude Z in the jet, a slice of the disk 
is seen under a solid angle d£l = Z dS /(r 2 + Z 2 ) 3/2 , where dS = 
dipirdr + Jr 2 /2) is the surface of the slice. Photons coming from 
this slice make an angle 6 S = arccos(Z/ Vr 2 + Z 2 ) with the axis 
(see figure 2). 

The luminosity of one face of the disk, for R in = 3 Rs and 
R out » Ri n , is given by the relation 

r Rout a Me 

Ldisk = crT disk (r) 2nr dr « —. (6) 
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Fig. 2: Disk radial and azimuthal splitting. A slice at (r, ip) e 
([Rin, Rout L [0, 2n]) is seen under a solid angle dfl from the jet at 
an altitude Z. 


Ttor (fit) 


D t (cos 6 t - a) 
IjurcKfitfi 


1/4 


Ldisk sin OJ 


0 


for cos 6 t e [-a : a] 
otherwise 


(7) 


with J(^) the distance between the slice center and the point-like 
source, 


d(fit) = A [(1 - a cos Of) 1 + a 2 sin 2 0,] / , (8) 

and tu the angle between the Z-axis and the emission direction 
from the disk: 


sin co = 


a sin# 

(1 + a 2 - 2acos6) l/2 


(9) 



Fig. 3: Dusty torus seen from an altitude Z in the jet under a 
solid angle dfl. The torus is sliced according to 6 t e \0 tmin ,6 tma J 
and tp € [0, 2n]. Each slice is illuminated by the disk and is in 
radiative equilibrium with a temperature T(fi t , ip) and emits as a 
blackbody (s = 1). 


3.3. Dusty torus 

The dusty torus is modeled by a torus shape structure whose 
major radius is D t and minor radius R t (see figure 3). 

In the same way as the disk, the torus is sliced in different 
parts that will radiate as blackbodies (so £ = 1 for the torus). 
Slices follow a linear discretization and are located with their 
coordinates (fi u ip) at the surface of the torus. 

Each slice of the torus is assumed to be in radiative equilib¬ 
rium with the luminosity received from the accretion disk. To 
simplify, we make the assumption that all the energy from the 
disk comes from its inner parts and that R in <tc ( D t - R t ) so that 
the source of energy is point-like when seen from the torus. With 

the parameter a = — < 1, the torus temperature is given by 


From an altitude Z in the jet, the torus is seen under a cer¬ 
tain solid angle which is delimited by 0 tmin and 0 tmax (see figure 
3). These values can be determined from geometrical considera¬ 
tions: 


6 tmin - arctan 


= arctan 




- arccos 


\ 

R t 


+ arccos 


\ 

Rt 

V z2+o ?, 


( 10 ) 


However, in the case of a continuum between the accretion disk 
and the dusty torus ( D t = R t + r out ) 9 6 tfnin will be chosen equal to 
0 . 


3.4. Broad line region 

The broad line region is modeled as an optically and geomet¬ 
rically thin shell of isotropically emitting clouds situated at a 
distance Rw r from the central black hole and extending up to an 
angle co max above the accretion disk plane (see figure 4). 



Fig. 4: The BLR, an optically and geometrically thin shell of 
isotropic clouds seen from an altitude Z in the jet under a solid 
angle d£l. The BLR is sliced according to co e [cL> max ,7 r/2] and 
ip e [0,2 n\. The BLR absorbs and re-emits part of the disk lumi¬ 
nosity. 

Tavecchio & Ghisellini (2008) showed that modeling the 
spectrum of the BLR with a blackbody spectrum at T = 10 5 /T 
provides a good approximation of the resulting inverse Comp¬ 
ton spectrum. We followed this idea using a temperature of 
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T b i r = 10 5 /T and an overall luminosity being a fraction a b i r of 
the disk luminosity. To achieve this, the BLR is modeled as a 
graybody at T b i r with an emissivity 

^blrLdisk /1 i \ 

&blr ~ ~ ~2 4 COS (O ma x- (H) 

^Klr^blr 

Thus the total luminosity of the BLR is given by 

Lblr — I &blri.U~Tblr^^ ~ ^blr^disk’ (12.) 

Jblr 

Like the torus, the BLR is divided linearly into slices along oj 
and cp. 

4. r eq in the jet 

4.1. Parameter values 

We have presented the description of the source modeling, and 
we can now choose the values for the different parameters. They 
are listed in Table 1. If not specified otherwise, these parameters 
are set for the rest of the study. Some characteristic values of the 
model are also derived in table 1. 


Parameter 

Symbol 

Value 

Black hole mass 

M 

5 x 10 8 M o 

BH accretion rate 

M 

1 M e dd 

Disk inner radius 

Rin 

3R S 

Disk outer radius 

Rout 

5 x 10 4 R S 

Disk emissivity 

Sdisk 

1 

Number of disk slices 

Nrdisk * A Pdisk 

18x3 

Torus center 

A 

10 5 R s 

Torus radius 

R t 

5 x 10 4 A 

Torus emissivity 


1 

Number of torus slices 

Nr tor X 

6x3 

BLR radius 

Rblr 

10 3 R s 

BLR angular opening 

COS CO m ax 

0.9 

BLR temperature 

Tbir 

10 5 7T 

BLR absorption 

Mblr 

0.1 

Number of BLR slices 

N r blr X N mr 

6x3 

Derived characteristic 

Symbol 

Value 

Schwarzschild radius 

Rs 

5.9 x 10 13 cm 

Disk temperature 

Tdisk 

[280 : 10 6 ] K 

Disk total luminosity 

Tdisk 

1.0 X 10 46 erg.s~ l 

Torus equilibrium 

Ttor 

[580 : 1345] K 

temperature 



Torus total luminosity 

Tfor 

1.4 X 10 45 erg.s~ l 

BLR total luminosity 

Tbir 

1.0 X 10 45 erg.s~ l 


Table 1: Parameters of the external sources. The values indicated 
in the right column are those used in this paper. 

With these parameters the external source spectra can be de¬ 
rived. An illustration of these spectra seen at an altitude Z = 
10 5 A in the jet is provided in figure 5. 

4.2. Resulting T eq 

We can now compute the resulting equilibrium bulk Lorentz fac¬ 
tor T eq all along the jet. This has been done for different set-ups 
of external sources (infinite accretion disk or finite accretion disk 
alone, finite disk + torus or finite disk + torus + BLR) and the 



log v [Hz] 

Fig. 5: Spectra of the three external soft photons sources seen at 
an altitude Z = 10 5 Rs on the jet axis in the sources frame. 



Fig. 6: Equilibrium bulk Lorentz factor resulting of the ex¬ 
ternal Compton emission for different external photon sources. 
The geometry is described in figure 1 with the following pa¬ 
rameters: finite and infinite accretion disk have an inner radius 
Ri n = 3 Rs ; the finite disk has an outer radius R out = 5 x 10 4 Rs ; 
Dtorus = 10 5 As, Rtorus = 5 x 10 4 R S , Rblr = 10 3 R S , and 
COS OJmax — 0.9. 


results are given in figure 6. In this plot, one can distinguish the 
effect of each external source. 

We start with the case of an infinite accretion disk ( R out = 
oo). For an emission zone in the jet, the inner parts of the accre¬ 
tion disk are always seen from below. As explained earlier, for 
a static source, this would lead to a Compton emission mainly 
directed toward the disk. This causes a transfer of momentum 
and a thrust forward on the plasma. However, as soon as the 
plasma accelerates, the photons coming from the outer part (for 
ji s < P~ l ) of the disk seem to travel backward in the comoving 
frame and produce a drag on the plasma. At every altitude, the 
equilibrium velocity is reached when the two effects balance. 
Analytical computation from Marcowith et al. (1995) showed 
that in the case of an infinite accretion disk, one should have 

Y eq = 1.16(|) 1/4 . With Rt = 3 R s , this gives T eq = 0.88 ( 
which is in agreement with our numerical results (see bold solid 
line compared to crosses in figure 6). 
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If we now consider a finite accretion disk ( R out = 5 x 10 4 Rs), 
we note the same behavior: at low altitudes the disk seems to be 
infinite seen from the axis. Once an altitude Z > R out is reached, 
the drag effect from the outer parts of the disk ceases and then 
the entire disk will imply a thrust on the bulk. As long as the 
acceleration is effective, Y eq will follow a law in Z/R out (Renaud 
& Henri (1998)). 

We can have the same reasoning concerning the effect of the 
dusty torus. As seen from the jet axis, the radiation from the 
torus comes at greater angles than the one from the accretion 
disk. Therefore, when the plasma accelerates, the torus radia¬ 
tion seems to come forward, which will tend to drag the flow. 
Nevertheless, in the lowest altitude the accretion disk radiation 
dominates and the resulting T eq is unchanged from the previous 
case. It is only from a certain altitude (Z « 10 3 Rs in our study) 
that the effect of the torus radiation starts to dominate and that 
the flow will actually slow down. Of course, the equilibrium ve¬ 
locity will never reach zero, as the radiation in the lab frame is 
never isotropic and always has a preferred direction upward. At 
one point (Z « R t in our study), most of the radiation from the 
torus moves forward in the comoving frame. This leads to a 
thrust on the flow and T eq increases again with the same acceler¬ 
ating slope as in the finite accretion disk case. 

The BLR photon field shows the same kind of effects with 
a deceleration regime inside the BLR (from 10 2 R$ to 10 3 Rs in 
our study), followed by an accelerating regime once the bulk 
leaves the BLR. At some point, the torus photon field becomes 
predominant and controls Y eq as explained previously. 

We note that the computation is done for a hot electron- 
positron plasma in the Thomson regime. For very hot plasmas, 
KN corrections will affect the rate of momentum transfer and the 
bulk Lorentz factor will differ in a way that is difficult to predict. 
Indeed, the Compton rocket is less efficient in the KN regime, 
but photons coming from larger incident angles are more likely 
to be in the KN regime and are precisely the ones dragging the 
flow. 

Moreover, far from the external photon sources, the relax¬ 
ation time to the equilibrium will become larger than the dy¬ 
namical time z/c. At this point, the acceleration will stop, lead¬ 
ing to an asymptotic value of the bulk Lorentz factor. The blob 
will then follow a ballistic motion. However, The point where 
this decoupling occurs depends on the absolute luminosity of 
the disk and the average energy of the plasma. A study of this 
phenomenon is presented in section 6. 

5. Influence of parameters on the equilibrium bulk 
Lorentz factor 

In this section we study the influence of the model parame¬ 
ters on T eq with the model composed of the accretion disk, the 
dusty torus at thermal equilibrium, and the BLR. If not other¬ 
wise stated, the parameters keep the values given in the previous 
section (table 1). 

5.1. Influence of the BLR on Y eq 

5.1.1. Influence of the BLR opening angle co max 

The influence of the BLR opening angle on T eq is shown in 
figure 7. The smaller oj max , the bigger the BLR, and the stronger 
the effect. 



Z/R s 


Fig. 7: Equilibrium bulk Lorentz factor for several BLR opening 
angles cos oj max . Values of the other geometrical parameters can 
be found in table 1. 

When (jj max increases, parts of the BLR at small co are sup¬ 
pressed. For the plasma inside the BLR (Z < Rbi r ), radiation 
from these parts moved backward in the comoving frame. The 
suppression of this radiation means less dragging effect and thus 
a higher Y eq for the flow inside the BLR. Thus, the differences 
between the different opening angles are important at Z < R b i r 

However, for the plasma outside the BLR (Z > R b i r ), radia¬ 
tion from these parts moved forward in the comoving frame. The 
suppression of this radiation means less thrust on the plasma, but 
the radiation from the parts of the BLR at greater (o, which seems 
to move forward in the comoving frame, is still present and so is 
the dragging effect. This is why the differences between the dif¬ 
ferent cases at Z > R b i r are not so important. Much of the thrust 
is given by the disk itself, even at altitudes close to the outer 
border of the BLR. 

5.1.2. Influence of the BLR absorption a bir 



Fig. 8: Equilibrium bulk Lorentz factor for several BLR absorp¬ 
tion a. Values of the other geometrical parameters can be found 
in table 1. 

The effect of the BLR absorption a b i r is shown in figure 8. 
The effect is similar to the opening angle effect. The stronger 
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the absorption, the stronger the emissivity and thus the stronger 
the luminosity from the BLR, and consequently, the stronger the 
drag effect. 

However, unlike in the opening angle (o max case, differences 
in the acceleration regime outside the BLR are noticeable. In¬ 
deed, with a reduction of the absorption, the drag and the thrust 
change, which was not the case previously. 

5.1.3. Influence of the BLR radius 

The effect of the BLR radius is shown in figure 9. The BLR 
radius is now going from Ru r - 10 3 Rs to Rbi r = 10 4 Rs. Other 
parameters are fixed to the values given in table 1. 



Fig. 9: Equilibrium bulk Lorentz factor for several BLR radius 
Rbi r . Values of the other geometrical parameters can be found in 
table 1. 

The effect of the source radius sizes is a bit different as it 
increases the amplitude of the drag or thrust, but also shifts the 
different regimes in altitude . We note that when the BLR ra¬ 
dius increases, its acceleration zone moves to higher Z. At one 
point (when Rbi r tends to D t - R t ), the radiation from the torus 
dominates and controls T eq . 

5.2. Influence of the torus on Y eq 

We study here the influence of varying R t , but always assume 
that D t = r out + R t . The results are shown in figure 10. 

The torus acts farther in the jet than does the BLR. Its radius 
effect is similar to the BLR’s: the bigger the radius, the stronger 
the drag, and the lower in the jet it occurs. With a smaller torus, 
higher T eq are reached closer to the black hole. Therefore, the 
emission from the lowest parts of the jet (at subparsec scale) 
will be strongly influenced by the torus size through the induced 
Doppler boosting. 


R blr /R S 



10 1 10 2 10 3 10 4 10 5 10 6 
Z/R s 

Fig. 10: Equilibrium bulk Lorentz factor for several torus sizes. 
R t changes as a free parameter as D t = r out + R t . The values of 
the other parameters can be found in table 1. 

shown that most of the emission is emitted within a characteristic 
emission cone of aperture angle « 1 /T. This led to the idea that 
the same object seen from a different angle will show a differ¬ 
ent broadband spectrum and led to the AGN unification scheme 
(Blandford & Rees 1978, Orr & Browne 1982 and Barthel 1989). 

For a given function T eq (Z), it is possible to compute the 
function of the equilibrium bulk Doppler factor S eq (Z,p 0 b s ), 
which depends on the altitude and on the observer viewing angle. 
Figure 11 shows the function 6 eq (Z,/j 0 b s ) in false colors, corre¬ 
sponding to the T eq (see also figure 11) computed in section 4 
with the source parameters given in table 1. Figure 12 shows the 
same 5 eq as a function of the altitude, but for four chosen obser¬ 
vation angles. 



5.3. Observation angle and Doppler factor 

The relativistic bulk Doppler factor is defined as 

Sb = 


1 


T& (1 fib Pobs) 


(13) 


/v 


with p 0 b s = cos i Q bs (see figure 1 for a definition of i 0 bs)- As — 

v- 3 

is a relativistic invariant (Rybicki & Lightman 1979), the spe¬ 
cific intensity in the lab frame is given by I v = /' S 3 b . It can be 


Fig. 11: Upper panel: Equilibrium bulk Doppler factor S eq in 
the jet plane altitude Z/R s vs observational angle i Q b s ( p 0 bs = 
cos i 0 bs)- The color scale is shown on the right. The green dashed 
line represents the corresponding J3 eq . Bottom panel: T eq as a 
function of the jet altitude. The geometry is described in section 
3 and the source parameters are listed in table 1. 

An observer situated at a constant p 0 b s sees the emission 
along the jet modulated by S eq (Z). It can be seen in figures 11 
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Fig. 12: Examples of equilibrium bulk Doppler factor as a 
function of the altitude for different observational angles i Q b s 
(fiobs = cos i 0 b s ). The green solid line represents the correspond¬ 
ing Y eq . The geometry is described in section 3 and the source 
parameters are listed in table 1. 


and 12 that 6 eq (Z ) shows several extrema at a constant p 0 b s . This 
means that certain zones of the jet are preferentially seen de¬ 
pending on the jet viewing angle i 0 b S . 

A few remarks can be made for peculiar values of 6b. 

a. Sb = 1 is an important value for the observer because it marks 
the limit between an increase and a decrease in the observed flux 
compared to the flux in the comoving frame. As Y eq (Z) varies 
in the observer frame, so does S eq , and the same observer can 
be in the emission cone of certain parts of the jet and out of the 
emission cone of other parts of the jet. With the bulk Lorentz 
factor computed in section 4, a value of S eq = 1 is possible for 
0.52 < ji () bs < 0.96. 


b. Extrema of S eq are found at Z verifying 


36b 1 3T eq \ 

^ =0o ^ r ,l^) =0 - 

We are then left with two possibilities: 



The solutions to this equation correspond to the T eq ex¬ 
trema. In our case, they are approximately at Z = 134 Rs, 
Z = 656R S ,Z= 1.73 x 10 3 R s , and Z = 3.61 x 10 4 fl s . 




Oh 

dY 


= 0 


It can be shown that this always happens for /3b = p 0 bs, or equiv¬ 
alently for Sb = Yb. These 6b(Z) extrema can be assimilated to 
a way in or out of the emission cone by the observer at p 0 b s . 
In our particular case, figure 11 shows that this is possible for 
min (fi eq ) = 0.45 < p Q b s < 1 = max(j3 eq ). Moreover, it can be 
confirmed for the cases p Q b s = 0.6 - 0.9 - 0.997 that the alti¬ 
tudes where J3 eq = p Q b s in figure 11 correspond to the altitudes 
where S eq reaches an extremum and where S eq = Y eq in figure 
12 . 


However, Y eq extrema can correspond to 6 eq minima or max¬ 
ima depending on p 0 b s . In the case J3 eq < ji 0 b s (blazar-type ob¬ 
jects), 6 eq maxima are correlated with Y eq maxima. In conse¬ 
quence, for these objects an observer will preferentially see jet 
zones where Y eq is at a maximum since the jet emission will 
be more boosted. On the contrary, for (3 eq > p Q b s (radio-galaxy- 
type objects), S eq maxima are correlated with Y eq minima. This 
means that an observer will not see the zones of the jet that have 
the highest speed, but - on the contrary- the jet emission will be 
dominated by the slowest zones. 

It is also interesting to note that there is a class of objects that 
will present very low modulation of the jet emission along Z. 

36 eq 

These objects are characterized by p Q b s ~ j3 eQ and thus —— ~ 0 

oZ 

almost everywhere in the jet. An example of this is shown by the 
case p 0 bs = 0.9 in figure 12, where S eq is almost constant from 
Z = 10 Rs to Z = 10 5 Rs . Of course other sources of variations 
can still produce an important variability for these sources. 

Similarly, different processes could dominate at different al¬ 
titudes, only due to Y eq evolution. In particular, external and syn¬ 
chrotron self-Compton emissions do not have the same beaming 
pattern for a given Doppler factor as shown by Dermer (1995). 
This author showed that synchrotron self-Compton follows a 
general beaming pattern ex S 3+a (with a the energy spectral in¬ 
dex of the radiation) whereas external Compton follows a beam¬ 
ing pattern oc 6 A b +2a because the Comptonized photon field is 
isotropic in the plasma rest frame in the SSC case whereas it 
depends on 6b in the external Compton case. 

Moreover, the computation in Dermer (1995) assumed an ex¬ 
ternal isotropic radiation and a pre-assumed bulk Lorentz factor. 
Both of these assumptions are no longer valid in our framework. 
This could have some consequences on the beaming statistics of 
these objects, but an exhaustive study of these effects could not 
be done without a complete modeling of the jet, which is not the 
purpose of this paper. 


6. Dependence of r# on the energetics 

As stated previously, in the two-flow paradigm the Compton 
rocket process finds its energy in the turbulence from the outer 
MHD jet through the relativistic particle emission. It is therefore 
understandable that the energetics of the particles will limit the 
influence of the Compton rocket effect on the actual value of Y. 
To compute the actual value of T, one can solve the following 
differential equation (Renaud & Henri (1998)), 


3Y b (Z,y e ) F' z 1 


OZ 


with F' : 


( 1 + ^) 


(14) 


4 nH' 


' = J* y'm e c 2 n' e (y')dy'. 


and p 

Here, it can be seen that the complete calculation of T^(Z, y e ) 
depends on the particle energy distribution. For the sake of sim¬ 
plicity, here we choose a Dirac distribution: 


n e (y) = Ne 6 (y- y e ). 

In this case, the following results are obtained: 


cr T 87T 


N e y 2 e H’ 


(15) 

(16) 
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and 

p' = N e y e m e c 2 . (17) 

For bulk Lorentz factors close to the equilibrium value, H' 
can be evaluated with a linear expansion: 

dH' 

H\T b ) * H\T eq ) + — (r* - r eq ). (18) 

By definition H'(T eq ) = 0, and it can be shown that 

< 19 > 

Peq L eq 

Finally, the following differential equation is solved: 


ar b (z,y e ) 

dZ 


-;A» (r '' (Z ' r ' )_r ‘ t<Z)) 


( 20 ) 


with /(Z, y e ) 


3m e C 3 Plq^lq 
8 ncr T y e H 


+ 


1 

3r 2 

eq 


As an example, we solve this equation and compute the ac¬ 
tual value of for several values of y e in two different cases: the 
accretion disk alone and the complete case seen section 4.2 with 
the accretion disk, the dusty torus, and the broad line region. In 
this study, we suppose y e constant along the jet. Of course, this 
is a simplistic assumption and a complex evolution of y e along 
the jet can be expected in a more complete modeling, but as we 
will see, it has very little influence on the evolution of in the 
lowest part of the jet. 

Here the Compton rocket appears as a restoring force on the 
plasma with a stiffness constant --. The relaxation length 

KZ,y e ) 

towards the equilibrium value, /(Z, y e ), is inversely proportional 
to y e (i.e., the higher y e is, the stronger the force is and the longer 
the plasma will actually follow T eq (Z)) , but the coupling is ef¬ 
fective in the lowest parts of the jet, no matter the value of y e . At 
some distance (of the order of /(Z, y e )) 9 the Compton rocket force 
slowly stops acting on the plasma. The bulk Lorentz factor then 
reaches a final Lorentz factor T^ = Too and follows a ballistic 
motion. It can be seen that in both cases in figure 13, high values 
of Too can be achieved with reasonable values of y e . The value of 
Too also depends strongly on the source geometry. Without the 
drag from the torus, it is easier for the jet to reach higher values 
of Too. Nevertheless, even with a strong torus, the computed val¬ 
ues of Too are entirely compatible with observed values by Lister 
et al. 2013 for the highest values of y e that are totally compatible 
with the observed high-energy emission. 




Fig. 13: Actual value of as a function of the altitude in the 
jet Z/Rs for several values of y e . Top: Standard accretion disk 
alone with an outer radius of Ae5R s . Bottom: Accretion disk, 
dusty torus, and broad line region with parameter values from 
table 1. 


the jet, displaying associated peaks in the emission that we call 
echoes. To show and study these effects, we will set up a very 
simple jet model where we compute the synchrotron radiation 
(SYN), the synchrotron-self Compton (SSC), and the external 
Compton (EC) radiations. The following results are just exam¬ 
ples of variations. A different model, or a complete modeling of 
the jet, would obviously show different results, but we can expect 
the variations to keep the same general features. 


7. Variation in the emission 

The aim of this section is not to create a realistic model of a jet 
nor to explain all the variability of a single object, but to illustrate 
what type of variations would be induced by a variation of T eq . 
We assume a jet composed of spherical emitting zones called 
blobs moving forward (see, e.g., Katarzynski et al. 2001, Boute- 
lier et al. 2008, or Hervet et al. 2015 for blob-in-jet models). 
The emergence of a blob at the base of the jet would then corre¬ 
spond to a flare. In our framework, this blob moves forward in 
the jet with an imposed bulk Lorentz factor Y eq . Because of the 
variation of Y eq inducing variations of S eq for an observer, the 
blob emission will show some interesting changes. Therefore, a 
single flare will induce a complex variability as it moves along 


7.1. Jet modeling 

To compute the emission, we first need to model the jet. The jet 
radius is fixed at a constant value Rj et = 5 x IQ 2 Rs . The magnetic 

field follows a power law B = Bo (^:) with Bo = 6.8 x 10 _2 G 
for the whole study. These values ensure that the same global 
evolution is kept between the magnetic energy density and the 
photon field energy density so we can compare the SCC emis¬ 
sion and the EC emission along the jet. 

We have chosen a pile-up distribution for the energy distribu¬ 
tion. It has the advantage of presenting one parameter less than a 
power-law distribution and concurs better with the idea of parti¬ 
cles accelerated through the second-order Fermi process. It can 
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be written as 

= ^3 eXP (-|)‘ (21) 

Particles have to be energetic enough to explain the high- 
energy emission of AGNs, but we also assumed the Thomson 
regime to compute Y eq . The Thomson condition can be expressed 
as ye s ( 1 - cos#*) «: 1 in the bulk frame. Using the sources pa¬ 
rameters given in table 1, the computation of es (1 - cos Os ) along 
the jet in the rest frame of a flow at T eq gives a maximum value 
of 10 -4 for the photons coming from the BLR, of 10 -5 for the 
photons coming from the disk, and of 10~ 6 for the photons com¬ 
ing from the torus. The pile-up distribution has a mean value 
< y >= 3 y and drops rapidly at higher y because of the expo¬ 
nential term. Setting a value of y - 10 5 allows us to stay in the 
Thomson regime almost everywhere along the jet and to obtain 
a Too > 5. Therefore, one can presume T = Y eq to compute the 
emission at every altitude in the jet. 

Pursuing what was said in the previous sections, each slice of 
the external emitting regions is a different source characterized 
by four numbers: p s = cos 0 S with 6 S the incoming angle and 
dQ the solid angle both described in figure 1 and figure 2; the 
temperature T ; and the emissivity s. These parameters seen in 
the bulk frame depend on ^ = T( 1 -J3bfi s ) and are given by 
(parameters in the bulk rest frame are denoted by a prime) 


II 

, 1 


d£Y = —d Q 

s l 

(22) 

, H~fo 

r 



£ = £ 



7.2. Computation of the emission along the jet in two energy 
bands 

The model being set, we were able to compute broadband emis¬ 
sion (including SYN, SSC, and EC) at every altitude along the 
jet for a flow at T^ = Y eq . From this, we computed the total emit¬ 
ted power at every altitude in two characteristic energy 

bands: 


- Infrared: [ImeV - leV] 

- y-rays: [20MeV - 300GeV] 


This allows us to study the evolution of the emission as a blob 
of particles moves forward in the jet. Nevertheless, the emission 
^ depends on T^ but also on the jet model. In order to decouple 
the variations of emission due to the jet model from the emission 
due to the variations of Y eq , we computed ( dP / dn) eq , be., the emis¬ 
sion with Yb = Y eq , and ( dP /do) 5 , he., the emission with Yb = 5 


0 dP/dfl) eq 

. This way, the function —— (Z) is only modulated by the 

(d P/ di2) 5 

variations of Y eq , excluding other sources of variability. 


of a single blob moving along the jet, the relation between Z and 
the observation time is given by 

J pZf ^ 

--(1 ~Peq A lobs) dZ- (23) 

Z„ PeqC X ’ 

An example of the dependance of the altitude on the ob¬ 
served time is given in figure 14 for several observation angles. 
The parameters to compute Y eq are the same as those used sec¬ 
tion 5.3. 



10 
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Fig. 14: Example of evolution of the observation time as a func¬ 
tion of the altitude in the jet for different observation angles. The 
solid line is the Y eq curve, the value of T is in the scale shown on 
the right. 


7.4. Emission variability as a function of model parameters 

In order to compare the evolution of the emission to the evolution 
of Y eq = along the jet, the study has been performed with the 
parameters given in table 1. The results are shown in figures 15 
and 16. 

The variations are simultaneous in infrared (IR) and in y rays 
and follow the variations of Y eq studied in section 5. The conclu¬ 
sions regarding the variations are quite similar to those on Y eq . 
The first echo lasting several hours is due to the acceleration of 
the flow followed by a deceleration due to the dragging effect 
from the BLR. The jet is then reaccelerated by the disk and the 
BLR before being dragged again by the torus, giving a second 
echo. 

However, as the flow moves more quickly, the time contrac¬ 
tion increases resulting in different variation timescale. These 
timescales depend on the sizes of the different sources of exter¬ 
nal emission. Because it is closer to the base of the jet, the BLR is 
responsible for the short timescale variations (from some hours 
to some days in our study). The torus, however, is responsible 
for variations at larger timescales (from several days to years). 


7.3. Time dependance of the emission in the case of a single 
blob moving in the jet 

In the case of a single blob traveling along the jet, we can convert 
the altitude into time, but owing to light time travel effects, the 
time between two events in the observer frame is different from 
the time between these two events in the lab frame. In the case 


7.4.1. Influence of the BLR parameters 

Figure 15 shows the influence of the BLR parameters on the time 
lag effects. The upper plot in figure 15 concerns the geometrical 
repartition of the BLR (through its opening angle oj max ) at con¬ 
stant total luminosity Lbi r = 0.1 . One can see that effects are 

more important with a larger covering factor because parts of the 
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(dP/dQ) eq 

(dP/dQ ) 5 


cos(a» max )=0.98 

cos(co max )=0.4 - 

cos (®max)= 0 - 9 . 

COS(co ma x)=0-1 . 

cos(co max )=0.7 . 




0 1 10 100 


Observation time (days) 


(dP/dQ) eq 

(dP/dQ ) 5 



«blr=° - 

a b |r=0.1 - 

ablr=°- 001 . 

a b |r=0.9 - 

a blr=°- 01 . 



0 1 10 100 


Observation time (days) 


(dP/dQ) eq 

(dP/dQ ) 5 


R b lr/R s = 1000.0 — 

R b ,r/R s = 7500.0 - 

R blr /R s = 2500.0 . 

R b , r /R s = 10000.0 - 

R b , r /R s = 5000.0 . 




Fig. 15: Evolution of the emission in the direction of the ob¬ 
server (cos i () bs = 0.997) as a function of observed time for a 
range of BLR parameters. The emission of the blob following 
r eq , (dP/d£l) eq is normalized by the emission of a blob at T = 5, 
(dF/dD) 5 . Values of all parameters are listed in table 1. The color 
scale indicates the dominant emission process. 


BLR closer to the jet axis have more influence on the Compton 
rocket. 

The influence of total luminosity of the BLR (given by abir) 
is more important (middle panel). For an ineffective BLR {au r - 
0), there is only one echo around day 4 imposed by the torus, 
but as the BLR becomes more important, two echoes appear, the 


first lasting several hours and the second a few days. We note 
that the greater the BLR, the more separated these two echoes 
are and the more peaked they are. Whereas the first echo always 
occurs around day 1, the second occurs later and later, from a 
few days to a hundred days. 

Last panel concern the BLR radius. The radius has almost the 
same effect as abir as it can delay the second echo. It is also worth 
noting that the first and second echoes are inversely important. 
As the BLR grows bigger, the first echo arises later and is more 
important (because T eq is). The second echo also arises later but 
because it is limited by the influence of the torus, its amplitude 
is diminished. 

Depending on the geometry and on the total luminosity of 
the BLR, we see that different behaviors in the time modulation 
of the emission are possible, which could lead to very different 
time variability in different objects. 

7.4.2. Influence of the torus parameters 

We can see the influence of the torus size in figure 16. The radius 
of the torus evolves while we keep the continuity between the 
accretion disk and the torus so D t = r out + R t . 



Fig. 16: Evolution of the emission in the direction of the ob¬ 
server (cos i a bs = 0.997) as a function of observed time for a 
range of torus parameters. The emission of the blob following 
r eq , (dP/dQ) eq is normalized by the emission of a blob at T = 5, 
(d/7dQ)5. Values of all parameters are given in table 1. The color 
scale indicates the dominant emission process. 

Here, we see again the two echoes imposed by the BLR and 
the drag from the torus after the flow crossed the BLR. As we can 
see, the second echo (which has timescale of at least ten days), 
is driven by the torus size. With a greater torus size, the echo 
arises sooner, but is more tamed. On the contrary, a smaller torus 
allows the flow to reach a larger velocity, implying a stronger 
echo here. 

At a certain point, as explained in section 2, the flow only ac¬ 
celerates, increasing the emission, but when the observer leaves 
the emission cone (which is highly dependent on the observation 
angle), the observed emission will decrease, giving the last echo. 
Here too, the torus size has a huge impact on the timescale of 
the variability. The smaller the torus, the sooner this echo arises 
(around day 100 here). When the torus size increases, this last 
echo gets delayed (up to day 2000 here), but its extent does not 
increase accordingly which makes it comparatively steeper. It is 
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also interesting to note that the maximum observed value of the 
emission does not depend on the torus size for this kind of echo. 

8. Conclusion 

The question of the acceleration of AGN jets is still a matter 
of discussion as we do not know the underlying processes or 
the precise speeds of the flows. The solution implied by the 
Compton rocket effect, viable in the two-flow paradigm, is 
elegant as it can naturally lead to relativistic speeds. In this 
work, we embrace this framework and study the influence of 
several external photon sources (the accretion disk, the dusty 
torus, and the broad line region) on the Compton rocket effect 
and on the induced bulk Lorentz factor. To do so, we carefully 
computed the resulting equilibrium bulk Lorentz factor, T eq , of 
a flow driven by the Compton rocket effect taking into account 
the anisotropy of the emission. With several external sources, 
T eq will show important changes along the jet, leading to ac¬ 
celeration and deceleration phases. We studied the influence of 
the external sources on these patterns and the induced Doppler 
factor as a function of the observation angle. We also showed 
that the emission of a flow following this T eq will experience 
correlated variations and that a single flare could be echoed 
several times. This could take part in the time-variation of these 
very variable objets and so we computed some examples of 
observed emission to illustrate our discussion. Even though 
these effects could not explain all the AGN variability alone, 
we find some interesting and non-trivial effects that could be 
part of the observed variability. This work could have some 
influence on the statistical study of AGNs and on their modeled 
variability. It will be followed by a more complete and more 
realistic model, applied to real objects to understand more 
precisely the role of the Compton rocket in the AGN variability. 
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